(19) 



J 



(12) 



(45) Date of publication and mention 
of the grant of the patent: 
06.12.2000 Bulletin 2000/49 

(21) Application number: 94111471.2 

(22) Date of filing: 22.07.1994 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets (11) EP 0 637 803 B1 

EUROPEAN PATENT SPECIFICATION 

(51) Intel. 7 : G06F 17/17, G06F 19/00 



(54) Method and device for adaptively estimating a transfer function of an unknown system 

Verfahren und Vorrichtung zur adaptiven Schatzung einer Obertragungsfunktion eines unbekannten 
Systems 

MSthode et dispositif d'estimation adaptatif d'une fonction de transfert d'un systeme inconnu 



(84) Designated Contracting States: 
DE FR 

(30) Priority: 27.07.1993 JP 18474293 

(43) Date of publication of application: 
08.02.1995 Bulletin 1995/06 

(73) Proprietor: 

NIPPON TELEGRAPH AND TELEPHONE 

CORPORATION 

Shinjuku-ku, Tokyo 163-19 (JP) 

(72) Inventors: 

• Tanaka, Masashi 
Tokorozawa-shi, Saitama (JP) 

• Kaneda, Yutaka 
Tokorozawa-shi, Saitama (JP) 

• Makino, Shoji 
Machida-shi, Tokyo (JP) 

• Haneda, Yoichi 
Nerima-ku, Tokyo (JP) 

• Kojima, Junjl 
Nerima-ku, Tokyo (JP) 



m 

CO 

o 

00 



CO 
CO 



GL 
LU 



(74) Representative: 

Hoffmann, Eckart, Dipl.-lng. et al 
Patentanwalt, 
Bahnhofstrasse 103 
82166 Grafelfing (DE) 

(56) References cited: 

• PROCEEDINGS OF ICASSP 89, IEEE NEW YORK 
USA, vol. 2, 23 - 26 May 1989, GLASGOW, UK, 
pages 1059-1062, XP002026916 T.FURUKAWA 
ET AL: 'The orthogonal projection algorithm for 
block adaptive signal processing" 

• IEICE TRANSACTIONS ON FUNDAMENTALS OF 
ELECTRONICS, COMMUNICATIONS AND 
COMPUTER SCIENCES, vol. E75-A, no. 11, 1 
November 1992, pages 1570-1578, XP000337105 
SUN W ET AL: "A NEW METHOD FOR 
PARAMETER AND INPUT ESTIMATION ON 
NONMINIMUM PHASE SYSTEMS" 

• PROCEEDINGS OF THE SOUTHEAST 
CONFERENCE (SOUTH EASTCON), 
BIRMINGHAM, ALABAMA, APR. 12 - 15, 1992, 
vol. 2, 12 April 1992, INSTITUTE OF 
ELECTRICAL AND ELECTRONICS ENGINEERS, 
pages 683-686, XP000338951 ZHOU P ET AL: 
"SPECTRAL ESTIMATION OF SEGMENTED 
SIGNAL USING DECONVOLUTION TECHNIQUE" 



Note: Within nine months from the publication of the mention of the grant of the European patent, any person may give 
notice to the European Patent Office of opposition to the European patent granted. Notice of opposition shall be filed in 
a written reasoned statement. It shall not be deemed to have been filed until the opposition fee has been paid. (Art. 
99(1) European Patent Convention). 

Printed by Xerox (UK) Business Services 
2.167 (HRS)/3.6 



EP 0 637 803 B1 

Description 

BACKGROUND QF THE INVENTION 

5 [0001] The present invention relates to a method for adaptivety estimating, with a projection algorithm, a transfer 
function of an unknown system and its output in an acoustic canceller, active noise control or the like and an estimating 
device using such a method. 

[0002] In the following description, time will be represented by a discrete time k. For example, the amplitude of a 
signal x at time k will be expressed by x(k). Fig. 1 is a diagram for explaining the estimation of the transfer function of 

10 an unknown system. Reference numeral 11 denotes a transfer function estimation part and 12 the unknown system, 
and reference character x(k) represents an input signal to the unknown system and y(k) an output signal therefrom. The 
transfer function h(k) of the unknown system is estimated using the input signal x(k) and the output signal y(k). Fig. 2 is 
a diagram for explaining an adaptive estimation of the transfer function. Reference numeral 21 denotes an estimated 
transfer function correcting vector calculation part, 22 an estimated transfer function correction part and 23 a convolu- 

15 tion part, these parts constituting an estimated signal generation part 20. The transfer function h(k) of the unknown sys- 
tem 12 is estimated as a transfer function h (k) of an FIR filter of a tap number L which forms the convolution part 23. 
More specifically, coefficients h ^k), .... h L (k) of the FIR filter are estimated. Let it be assumed that the "transfer func- 
tion" and the "FIR filter coefficient" will hereinafter be construed as the same. For the sake of brevity, the filter coefficient 
is represented as an estimated transfer function vector h (k) defined by the following equation. 

20 

h(k) = [h 1 (k),h 2 (k) h L (k)] T (1) 

where T represents a transpose. 

[0003] In Fig. 2, the input signal x(k) to the unknown system 12 is fed to the convolution part 23 and a calculation 
25 is performed to obtain the estimated transfer function vector h (k) that minimizes an expected value which is the square 
of an error signal e(k) available from a subtracter 24 which detects a difference between an output y (k) from the con- 
volution part 23 given by the following equation (2) and the output y(k) from the unknown system 12. 

y(k) = h(k) T x(k) (2) 

30 

x(k) = [x(k),x(k-1) x(k-L+1)] T (3) 

where y (k) is an estimated value of the output from the unknown system, which is close to the value of the output y(k) 
when the estimated transfer function h (k) is close to an unknown characteristic. 

35 [0004] In practice, the transfer function of an unknown system often varies with time as in the case where the trans- 
fer function of an acoustic path varies with movement of audiences or objects in a sound field such as a conference hall 
or theater. For this reason, the estimated transfer function h (k) of the unknown system also needs to be adaptively cor- 
rected in accordance with a change in the acoustic path of the unknown system. The estimated transfer function cor- 
recting vector calculation part 21 calculates an estimated transfer function correcting vector 8h (k) on the basis of the 

40 error signal e(k) and the input signal x(k) to the unknown system 12. The estimated transfer function correction part 22 
corrects the estimated transfer function by adding the estimated transfer function correcting vector 6h (k) to the esti- 
mated transfer function vector h (k) as expressed by the following equation. 

h(k+1) = h(k) + u5h(k) (4) 

45 

where u is called a step size, which is a preselected quantity for controlling the range of each correction and is handled 
as a time-invariant constant. In the following description, the estimated transfer function correcting vector 5h (k) is cal- 
culated for u = 1 and, if necessary, it is multiplied by a desired step size u.. In some applications the value of the step 
size u is caused to vary with time, but in such a case, too, the following description is applicable. The corrected esti- 
50 mated transfer function vector is transferred to the convolution part 23. The above is a transfer function estimating oper- 
ation at time k and the same operation is repeated after time k + 1 as well. 

[0005] T. Furukawa et al: The Orthogonal Projection Algorithm for Block Adaptive Signal Processing" in Proceed- 
ings of ICASSP 98, IEEE New York USA, Vol. 2, 23-26 May 1989, Glasgow UK, pp. 1059-1062, disclose estimation of 
a transfer function of an unknown system by using orthogonal projection based on an Affine Projection Algorithm. The 
55 aim is to reduce the number of multiplications, wherein a filter coefficient vector is updated for each block of input vec- 
tors (or every r sample times) to reduce computational requirement. Since the coefficient vector is updated for each 
block L, the number of update times is reduced, resulting in decrease in convergence speed of the adaptive filter. The 
filter coefficient vector is adjusted according to an equation numbered 1 1 in the document and including the Moore-Pen- 
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rose type generalized inverse matrix which still involves a high computational requirement. The projection algorithm 
proposed by T. Furukawa et al. is developed from the Affine Projection Algorithm proposed by Ozeki and Umeda and 
mentioned below. 

[0006] The estimated transfer function correcting method described above in respect of Fig. 2 is known as an adap- 
5 tive algorithm. While an LMS (Least Mean Square) algorithm and an NLMS (Normalized LMS) are also well-known as 
adaptive algorithms, a description will be given of a projection algorithm proposed in a literature [Ozeki and Umeda: "An 
Adaptive Filtering Algorithm Using an Orthogonal Projection to an Affine Subspace and Its Properties", Journal of Insti- 
tute of Electronics, Information and Communication Engineers of Japan (A), J67-App. 126-132, (1984-2).] This docu- 
ment discloses a method of adaptively estimating thetransfer function of an unknown system in accordance with the first 
10 portion of claims 1, 2 and 3. 

[0007] The projection algorithm requires a larger number of operations than does the NLMS but has an excellent 
adaptability to a speech signal input. With the projection algorithm, as referred to above, the vector 6h (k) is determined 
by Eq. (4) for u. = 1 so that simultaneous equations composed of the following p equations are satisfied. 

15 y(k)- = x(k) T (h(k) + 8h(k>) 

y(k-l) = x(k-l)T(h(k> ♦ &(k)) 

y(k-p+l) = x(k-p+l) T (h(k) + Sh(k) ) (5) 

Eq. (5) indicates that the vector 8h (k) is determined so that the estimated transfer function h (k+1) updated at time k 

provides the same values y(k), y(k-1) y(k-p+1) as the outputs from the unknown system, respectively, for p input 

25 vectors x(k), x(k-1) x(k-p+1) prior to time k. By this, it is expected that the characteristic of the estimated transfer 

function h (k) will approach the characteristic of the unknown system as the adaptive updating of the estimated transfer 
function is repeated using the vector 5h (k). 

[0008] In the above, p is a quantity commonly called a projection order. As the projection order p increases, the 
adaptability of the projection algorithm increases but the computational complexity also increases. The conventional 
30 NLMS method corresponds to the case of p = 1 . 

[0009] Now, transposing the first equation in Eq. (5), we have 

x(k) T 5h(k) = y(k) - x(k) T h(k) = e(k) (6) 



35 Furthermore, transposing the second equation in Eq. (5) and using an equation obtainable by setting k in Eqs. (4) and 
(6) to k-1, we have 



40 



x(k-1) T 5h(k) = y(k-1)-x(k-1) T h(k) 

= y(k-1) - x(k-1) f h(k-1) + u5h(k-1)) 

= y(k-1) - x(k-1) T h(k-1) + ux(k-1) f 5h(k-1) 

= e(k-1)-ue(k-1) 

= d-u)e(k-1) 



(7) 



45 Thereafter, the following relationship similarly holds. 

x(k-i) T 8h(k) = (1-u) l e(k-i) (8) 
Based on this, Eq. (5) can be rewritten as the following system of simultaneous equations. 

50 

X p (k)5h(k) = e(k) (9) 

where X p (k) is a matrix with p rows and L columns and e(k) is a vector of the p rows; they are defined by the following 
equations. 
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X p OO = 



x(k)* 
x (k-l) T 



(10) 



10 



15 



• (k) 



e(k) 
(l-Me(k-l) 



(i-n) p_1 e(k-p+i)_ 



(11) 



20 The vector e(k) will hereinafter be referred to as an error vector. Now, since p is usually smaller than L, Eq. (9) is an 
under-determined simultaneous equation or indeterminate equation for the vector Sh (k) and the minimum norm solu- 
tion of the vector 5h (k) is given by the following equation. 



25 



5h(k) = X D (k) T (X D (k)X D (k) T )- 1 e(k) 



= [x(k)x(k-1),. 



•pw-pv 
,x(k-p+1)]g(k) 



(12) 



where 



30 



35 



g(k) = R p (k)- 1 e(k) 
R p (k) = X p (k)X p (k) T 



(13) 
(14) 



R p is a matrix with p rows and p columns, which will hereinafter be referred to as a p-order covariance matrix or auto- 
correlation matrix, and g(k) a pre-filter coefficient vector. Letting elements of the pre-filter coefficient vector g(k) be rep- 
resented by g^k), g 2 (k), .... g p (k), the estimated transfer function correcting vector 8h (k) can be expressed on the basis 
of Eq. (12) as follows: 



40 



5h (k) = £9i(k)x(k-i + 1) 



(15) 



45 



50 



55 



Thus, when the projection algorithm is used, in accordance with the present invention the estimated transfer function 
correcting vector calculation part 21 in Fig. 2 has such a construction as depicted in Fig. 3. Reference numeral 31 
denotes a pre-filter coefficient vector calculation part, which uses the input signal x(k) and the error signal e(k) to cal- 
culate the pre-filter coefficient g(k) on the basis of Eq. (13). Reference numeral 32 denotes a pre-filtering part, which 
performs the pre-filtering operation expressed by Eq. (15) to synthesize the estimated transfer correcting vector 8h (k) 
by use of the pre-filter coefficient g(k) that is transferred from the pre-filter coefficient vector calculation part 31. 
[0010] Now, a description will be given of the computational complexity of the projection scheme described above. 
The computational complexity mentioned herein is the number of multiplication-addition (or addition) operations neces- 
sary for estimating operations per unit discrete time. The computational complexity of Eq. (2) in the convolution part 23 
of the tap number L in Fig. 2 is L. The computational complexity of Eq. (13) in the pre-filter coefficient vector calculation 
part 31 of the estimated transfer function correcting vector calculation part 21 is about p 3 /6 when using the Choleski 
method which is a typical computation method. The computational complexity of Eq. (15) in the pre-filtering part 32 is 
(p-1)L The computational complexity of Eq. (4) in the estimated transfer function correction part 22 is L. Thus, the entire 
computational complexity NC of the projection scheme is given by the following equation. 
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NC = L + p 3 /6 + (p-1)L + L (16) 

[0011] On the other hand, the computational complexity of the NLMS scheme or LMS algorithm is about 2L For 
example, when L = 500 and p = 20 (a typical value in the case of an acoustic echo canceler), the number of operations 
5 involved in the NLMS scheme is 1000, whereas the projection scheme requires as many as about 12000 operations on 
the basis of Eq. (16). The computational complexity p 3 /6 of Eq. (13), in particular, abruptly increases as the projection 
order p increases. Thus, the projection scheme has excellent convergence characteristics as compared with the NLMS 
scheme but poses the problem of increased computational complexity. 

10 SUMMARY OF THE INVENTION 

[0012] It is therefore an object of the present invention to provide an adaptive transfer function estimating method 
which permits reduction of the computational complexity involved and an estimating device using such a method. 
[0013] This object is achieved with a method as claimed in claim 1 or 2 and a device as claimed in claim 7 or 8. 
15 Preferred embodiments are subject-matter of the depnedent claims. 

[0014] To attain the above object, a first embodiment of the present invention reduces the computational complexity 
of the pre-filter coefficient vector calculation part 31 while a second embodiment and/or reduces the computational 
complexity of the pre-filtering part 32. 

20 PRIEF DESCRIPTION OF THE DRAWINGS 
[0015] 

Fig. 1 is a block diagram illustrating an ordinary construction for the estimation of a transfer function; 
25 Fig. 2 is a block diagram of a transfer function estimation part 1 1 in Fig. 1 ; 

Fig. 3 is a block diagram of an estimated transfer function correcting vector calculation part 21 embodying the 
present invention in Fig. 2; 

Fig. 4 is a block diagram of a pre-filter coefficient vector calculation part 31 according to the first aspect of the inven- 
tion in Fig. 3; 

30 Fig. 5 is a block diagram showing an estimated signal generating process according to the second aspect of the 
invention; 

Fig. 6 is a block diagram showing an application of the present invention to measurement of the transfer function of 
a loudspeaker; 

Fig. 7 is a block diagram showing an application of the present invention to an echo canceller; 
35 Fig. 8 is a block diagram showing an application of the present invention to noise control; 

Fig. 9 is a graph showing echo cancelling characteristics of an echo canceller embodying the present invention; 
and 

Fig. 10 is a graph showing, in comparison with a prior art example, the relationship between the number of opera- 
tions for the estimation of the transfer function and the projection order. 

40 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0016] A description will be given first of a method for reducing the computational complexity of the pre-filter coef- 
ficient vector calculation part 31 . This method is characterized in that the pre-filter coefficient vector g(k) is obtained by 
45 a recursion formula utilizing a vector g(k-1) at one unit time before. Now, it will be demonstrated how the vector g(k) can 
be obtained from the immediately preceding vector g(k-1) on a recursive basis. The recursion formula is derived 
through utilization of the fact that elements of the error vector e(k) of Eq. (11) shift with time while being multiplied by 1- 
H and the property of the inverse of the covariance matrix of the input signal given by Eq. (14). 
[0017] As is the case with Eq. (14), a covariance matrix Rp.i(k) of the p-1 order is defined by the following equation. 

so 



55 
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Vi (k) s x p-i (k> Vi (k) 



10 



Mk) 
(k-1)' 



x<k-p+2) J 



[x (k),x (k-1) 



x (k-p+2)] 
(17) 



J5 



20 



25 



[0018] It is known from a literature, S. Haykin, "Adaptive Filter Theory," 2nd edition, Prentice-Hall, 1991, pp. 577- 
578 that the inverse of the covariance matrix Rp(k) and the inverse of the covariance matrix Rp.i(k) bear such a rela- 
tionship as shown below. 



Rp(k) 



-1 



Rp-I<l<-1> 



-1 



a(k)a(k) 
Ftk) 



(18) 
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35 



Rp(k-l)" 1 



-1 



. b(k-1)b(k-1 )T (ig) 
B(k-l) 
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where a(k) is a forward linear prediction coefficient vector (p-order) which satisfies a normal equation 
R p (k)a(k) = [F(k), 0, .... 0] T and its first element is a 1; F(k) is the minimum value of the sum of forward a posteriori pre- 
diction-error squares; b(k-1) is a backward linear prediction coefficient vector (p-order) which satisfies a normal equa- 
tion R p (k-1)b(k-1) = [0 0, B(k-1)] T and its last element is a 1; and B(k-1) is the minimum value of the sum of 

backward a posteriori prediction-error squares. The linear prediction coefficients a(k), b(k-1) and the minimum values 
of the sums of forward and backward a posteriori prediction-error squares F(k) and B(k-1) could be calculated with less 
computational complexity through use of an FTF (Fast Transversal Filters) algorithm, for instance, as disclosed in J. M. 
Cioffi and T. Kailath, "Windowed fast transversal adaptive filter algorithms with normalization," IEEE Trans. Acoust, 
Speech Signal Processing, vol. ASSP-33, no. 3, pp. 607-625. The pre-filter coefficient vector g(k) of Eq. (13) is shown 
again below. 



g(k) = R p (k)- 1 e(k) 



(20) 



A pre-filter deriving coefficient vector f(k) of p-1 order is defined, corresponding to the vector g(k), by the following equa- 
tion. 



f(k) = R M (k)" , e M (k) 



(21) 



6 



where 



EP 0 637 803 B1 



P-l 



(JO 



10 



e(k) 
(l-Ji)e(k-l) 



(l-ji) p ~ 2 e(k-p+2)_ 



(22) 



Since all elements of a vector (1-n)e p _ 1 (k-1) of p-1 order, obtained by substituting (k-1) for k in Eq. (22) and multiplying 
both of its left and right sides by constitute elements of the p-order vector e(k) of Eq. (1 1) except its first element 
15 e(k), the relationship of the following equation holds. 



20 



(k) 



!(k) 



(1.(1). (k-1) 
p-l 



(23) 



25 Furthermore, since all the elements of the p-1 order vector of Eq. (22) constitute all elements of the p-order vector of 
Eq. (11) except its last element, the relationship between the vectors e(k) and e^k) can also be expressed as follows: 



30 



e (k) = 



' -p-1 00 

(1- p. ^"^(k-p+l) 



(24) 



35 



40 



Multiplying the respective terms on both sides of Eq. (18) by e(k) from the right, we obtain the following equation from 
Eqs. (20), (21) and (23). 



*(k) = (1- 11) 



\° 1 



+ a<k) e(k) a(k) (25) 
F(k) 



45 



Moreover, multiplying both sides of Eq. (19) by e(k-1) from the right, we obtain the following equation from Eqs. (20), 
(21) and (24). 



so 



™-[ ic r] 



+ b(k-D •<k-i) b(k _ 1) 

B(k-l) 



(26) 



Transposing the right side to the left side, we have 
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P <k " 1 'l - . - b ,27, 

L 0 J B(k-l) 



In this way, the pre-filter coefficient vector g(k) is calculated on the basis of the value of the pre-filter deriving coefficient 
w vector f(k-1) as expressed by Eq. (25), and the vector f(k-1) is calculated from the vector g(k-1) as expressed by Eq. 
(27). That is, Eqs. (25) and (27) are recursion formulae for the vector g(k). 

[0019] On the right side of Eq. (25), to obtain the first term (1-u)[ ] needs p-1 multiplications to obtain e(k) needs p- 
1 multiplications and to obtain the second term a(k) T e(k) needs p multipty-add operations, one division for division by 
F(k) and p multiply-add operations for multiplication by a(k) and addition to the first term on the right side. That is, a total 

is of 4p operations are needed. Similarly, approximately 2p operations are required to obtain Eq. (27). Accordingly, the 
number of operations of Eqs. (25) and (27) is around 6p, which is a computational complexity proportional to p. On the 
other hand, it is disclosed in the afore-mentioned literature by J. M. Cioffi and T. Kailath that the linear prediction coef- 
ficient vectors a(k), b(k-1) and the minimum values of the sum of a posteriori prediction-error squares F(k), B(k-1), 
which are needed in Eqs. (25) and (27), can be calculated by a linear prediction scheme with a computational complex- 

20 ity of about 10p. Thus, it is evident that the recursion formulae of Eqs. (25) and (27) permits calculation of the vector 
g(k) with the computational complexity proportional to the projection order p. 

[0020] According to the present invention based on the above-described principles, the matters (A), (B) and (C) 
listed below are effective in reducing the computational complexity and in increasing the computational stability. 



25 (A) Eq. (27) is an equation for the p-order vector and its left side is zero with respect to a p-th element (the least 
significant element) of the vector. Since a p-th element of the vector b(k-1) is always a 1, the following equation 
holds using a p-th element g p (k-1) of the vector g(k-1). 

1 " k - 1 ' T " k - 11 . g.k-1, US, 
B(k-l) P 



35 Therefore, the following equation may be calculated in place of Eq. (27). 



[f (k-l)~| 



4 o | n | = (k-1) -g (k-1) b (k-1) (27a) 



In this case, however, taking computational errors into account, it may sometimes be advantageous to obtain both 
45 9p(k-1) and the left side of Eq. (28) and average them. 

(B) When denominators F(k) and B(k-1) of the second terms on the right sides of Eqs. (25) and (27) are small, oper- ' r 
ations are unstable. The computational instability could be reduced by adding non-negative 8 F (k) and 6 B (k-1) to the 
denominators as shown in the following equations. 
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[_f (k-l)J 



a (k)T « (k) a (k) (29, 



b(k l) T e(k-l) ■ (30) 



1 1 B(k-l) *5 (k-1) 

15 In practice, the values ^(k) and 6 B (k-1) may be set to desired values about 40 dB smaller than the average power 
of the input signal x(k) (that is. about 1/10000 the average power) or may also be varied with time k in accordance 
with power variations. 

(C) In linear prediction analysis for obtaining the linear prediction coefficient vectors a(k), b(k-1) and the minimum 
values of the sum of prediction-error squares F(k), B(k-1) which are needed in Eqs. (25) and (27), the analysis 

20 frame of the input signal x(k) differs from that in an ordinary case. In concrete terms, the analysis frame in an ordi- 
nary linear prediction analysis ranges from time 0 to the current time, and when time k-1 is updated to k, x(k) needs 
only to be added to the analysis frame. In the projection algorithm, the analysis frame ranges from x(k) to 
x(k-L-p+2) in Eq. (5); so that when time k-1 is updated to k, it is necessary not only to add x(k) to the analysis frame 
but also to remove x(k-L-p+1) . On this account, the linear prediction analysis in the projection algorithm requires 

25 computational complexity twice that needed for ordinary linear prediction analysis. However, when the statistical 
property of the input signal does not change or when it is expected that its change is slow, the result of the linear 
prediction analysis does not largely depend on the analysis frame; hence, it is possible to use the ordinary linear 
prediction analysis in which x(k) needs only to be added to the analysis frame at time k - this permits reduction of 
the computational complexity. 

30 

[0021] Next, a description will be given of an embodiment of the adaptive transfer function estimating method of the 
present invention based on the above-described theoretical discussions, reference being made to Fig. 2 because the 
overall construction of its functional block is the same as shown in Fig. 2. Since the above-described theory underlying 
the present invention is related to the reduction of computational complexity in the estimated transfer function correcting 
35 vector calculation part 21 in Fig. 2, in particular, in the pre-filter coefficient vector calculation part 31 and the pre-filtering 
part 32 depicted in Fig. 3, these parts will hereinbelow be described in detail. 

[0022] Fig. 4 illustrates in block form the pre-filter coefficient vector calculation part 31 based on the discussion 
above. Reference numeral 41 denotes a linear prediction part, 42 a pre-filter deriving coefficient vector correcting part, 
43 a pre-filter coefficient vector correcting part, and 44 an error vector generating part. The linear prediction part 41 cal- 

40 culates the forward linear prediction coefficient vector a(k) which satisfies a normal equation 
R (k)a(k) = [F(k), 0 0] T and the sum of forward a posteriori prediction-error squares F(k) when the prediction coef- 
ficient a(k), the backward linear prediction coefficient vector b(k) which satisfies a normal equation 
R (k-1)b(k-1) = [0 0, B(k-1)] T and the sum of backward a posteriori prediction-error squares B(k) when the predic- 
tion coefficient b(k) is used. These values can be calculated by methods disclosed in the afore-mentioned literature by 

45 J. M. Cioffi et al. The error signal vector generating part 44 stores p error signals e(k), e(k-1). .... e(k-p+1) and consti- 
tutes the error signal vector e(k) of Eq. (1 1). In the pre-filter deriving coefficient vector correcting part 42, the pre-filter 
coefficient vector g(k-1), the backward linear prediction coefficient vector b(k-1), the error signal vector e(k-1) and the 
minimum value of the sum of backward a posteriori prediction-error squares B(k-1) are used to compute the pre-filter 
deriving coefficient vector f(k-1) on the basis of Eq. (27) or (30). In the pre-filter coefficient vector correcting part 43, the 

so pre-filter deriving vector f(k-1), the error signal vector e(k), the forward linear prediction coefficient vector a(k) and the 
minimum value of the sum of forward a posteriori prediction-error squares F(k) are used to compute, by Eq. (25) or (29), 
the vector g(k) that satisfies Eq. (13). 

[0023] With the methods mentioned above, the computational complexity involved in the pre-filter coefficient vector 
calculation part 31 can substantially be reduced from p 3 /6 to 15p. Since the computational complexity (p-1)L in the pre- 
ss filtering part 32 remains unsolved, a large number of operations are still needed when L is large. 

[0024] Next, a description will be given of a solution to the above-noted problem by storing approximate values of 

the estimated transfer function vector h (k+1) and averaging pre-filtering coefficients. 

[0025] At first, substituting k-1 for k in Eq. (4), the estimated transfer function h (k) is expressed as follows: 
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h(k) = h(k-1) + u5h(k-1) 



(31) 



Substituting this in Eq. (4), the estimated transfer function h (k+1) is expressed as follows: 



h(k+1) = h(k-1) + iiSh(k) + fi8h(k-1) 



(32) 



Setting k-2, k-3, ... for k in Eq. (4) and repeating the above substitution, the estimated transfer function h (k+1) is 
expressed by 



In this instance, the estimated initial value h (0) is set to 0 and this equation reveals that the estimated transfer function 
h (k+1) is a summation of correcting vectors uSh (k), u8h (k-1), .... |i8h (0) from time 0 (the transfer function estimation 
starting time) to the current time k. 

[0026] The correcting vector uSh (k) is expressed by Eq. (15). Setting k-1, k-2, for k in Eq. (15) and substituting 
Eq. (15) in Eq. (33), we get 



h(k+1) = u[{ gi (k)x(k)+g 2 (k)x(k-1)+ ... + g p (k)x(k-p+1)} + { gi (k-1)x(k-1)+g 2 (k-1)x(k-2) 

+ ... +g p (k-1)x(k-p)} + ... + { gi (0)x(0)+g 2 (0)x(-1) + ... + g p (0)x(-p+1)}] 
= u[ gi (k)x(k) +{g 2 (k)+ gi (k-1)}x(k-1) +{g 3 (k)+g 2 (k-1)+ gi (k-2)}x(k-2) + ... 

+ {gp.i(k)+g p _ 2 (k-1)+ ... + g1(k-p+2)}x(k-p+2) + {g p (k)+g ^-1)+ ... + g 2 (k-p+2) 
+ g 1 (k-p+1)}x(k-p+1) + {gpfk-D+g^) + ... + g 2 (k-p+1)+ gi (k-p)}x(k-p) 
+ {g p (k-2)+g p _ 1 (k-3)+ ... + g 2 (k-p) + g ^k-p-IJMk-p-l) ♦ 



[0027] From this equation the following facts are expected. Firstly, the pre-filter coefficient gj(k) is calculated at 
every time k in the pre-filter coefficient calculation part 31 in Fig. 3 and provided to the pre-filtering part 32; in this case, 
it is expected that the computational complexity would be reduced by smoothing (or averaging) the pre-filter coefficient 
gi(k). Secondly, the term + g p (k-1)+ ... and the subsequent terms in Eq. (34) do not involve the pre-filter coefficient gtfk) 
which is fed at time k, and hence this portion does not change after time k. By storing these terms as approximate val- 
ues of the vector h (k+1), it is expected that the computational complexity would be reduced accordingly, because no 
calculations are needed for these terms after time k. 

[0028] Next, the above will be expressed by mathematical formulae. The smoothing of the pre-filter coefficient takes 
place for each corresponding input vector x(k-i). From Eq. (34), the smoothing corresponding to the input vector x(k-1), 
for example, is g 2 (k)+g n (k-1) , and the smoothing corresponding to the input vector x(k-2) is g 3 (k)+g 2 (k-1)+g j(k-2) . 
Letting the result of averaging (a smoothing coefficient) corresponding to the input vector x(k-i+1) , inclusive of the con- 
stant n, be represented by Sj(k), it is expressed by 



h(k+1) = ii8h(k) + n5h(k-1) + ... + u6h(0) 



(33) 



i-1 



Sj(k)=uXg H (k-j) forlsisp 



(35) 



j=o 



i-1 



Sj(k)=u]Tg 4 (k-j) forp<i 



(36) 



Eq. (35) is expressed as follows: 
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M 

Si(k) = ^g 4 (k-j) + ngi(k) (37) 

,2 jsl 
= ^Z9H-i( k + 1 ) + Wi( k ) 

= s M (k-1) + u gi (k) for 2 ^ i S p 
= u 9j (k) fori = 1 



[0029] On the other hand, letting the approximate value of the estimated transfer function h (k+1) be represented 
by z(k+1), it is expressed by 

15 

z(k+D = fapM+gp^k-l)* ... +g 2 (k-p+2)+g ^k-p+IJMk-p+l) 

+{g p (k-1)+g M (k-2)+ ... +g 2 (k-p+1)+g ^k-p^k-p) (38) 
+{g p (k-2)+g p _ 1 (k-3)+ ... +g 2 (k-p)+ gi (k-p-1)}x(k.p-1) + ... 

20 

From Eqs. (34), (35) and (38) we have the estimated transfer function h (k+1) expressed as follows: 

p-1 

25 h (k+1) = z(k+1) + £ Sj(k) x (k-i+1) (39) 

Moreover, the following relationship holds between the approximate values z(k+1) and z(k). 

30 

z(k+1) = {g p (k)+g H (k-1)+ ... +g 2 (k-p+2) +g ^k-p+IJMk-p+D+zM (40) 
= s p (k)x(k-p+1)+z(k) 

35 [0030] The estimated value y (k) of the output from the unknown system is expressed, from Eqs. (2) and (39), as 
follows: 



y(k) = x(k) T h(k) (41) 

40 

= x(k) T {z(k) + £ Sj (k-1)x(k-i» 

i=1 
p-1 

= x(k) T z(k) + Xsj(k.1)x(k) T x(k-i) 

45 i=1 

= x(k) T z(k) + s p . 1 (k-1) T r p . l (k) 

so In the above. Sp.^k-1) is a smoothing coefficient vector and r^k) is a correlation vector, which are defined by the fol- 
lowing equations. 

Sp.^k-1) = [s^k-1), s 2 (k-1) s^fk-l)] 7 (42) 

55 r^k) = (x(k) T x(k-1), x(k) T x(k-2) x(k) T x(k-p-1)] T (43) 

Since the vector x(k), defined by Eq. (3), is given by 
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x(k) = [x(k), x(k-1) x(k-L+1)] T (44) 

the following relationship holds 

5 x(k) T x(k-i) = x(k-1) T x(k-M) -x(k.L)x(k-L-i)+x(k)x(k-i) (45) 

where i = 1, 2 p-1 

and the following equation holds 

10 r M (k) = r M (k.1).x(k.L)x M (k-L) +x(k)x p . 1 (k) (46) 

where 

Xp^k) = [x(k-1), x(k-2) x(k-p+1)] T (47) 

15 

[0031] Next, a description will be given, with reference to Fig. 5, of the storage of the approximate value z(k) of the 
estimated transfer function h (k) and the transfer function estimation procedure which is followed when the pre-filter 
coefficients are smoothed. At time k, a correlation calculation part 52, which is supplied with the input signal x(k), cal- 
culates the correlation vector r^k) by Eq. (46), using the input signal x(k), previous input values x(k-1) x(k-L) and 

20 the correlation vector r^k-l) at the immediately preceding time. 

[0032] Then, an inner product s ^(k-1) r ^(k) of the correlation vector r^k) and the smoothing coefficient vec- 
tor s p_ t (k- 1 ) of the pre-filter coefficients is calculated in an inner product calculation part 53. A convolution part 54 per- 
forms a convolution x(k) T z(k) of the stored approximate value z(k) of the estimated transfer function and the input signal. 
The results of the inner product calculation and the convolution are added together in an addition part 57 to synthesize 

25 the estimated value y (k) of the unknown system output. These operations correspond to the operation of Eq. (41). 
[0033] Following this, the error e(k) between the unknown system output y(k) and the estimated value y (k) is 
obtained by the subtracter 24 shown in Fig. 2 and the pre-filter coefficient g ( (k) is calculated in the pre-filter coefficient 
vector calculation part 31 shown in Fig. 4. 

[0034] After this, the calculated pre-filter coefficients are sent to a pre-filter coefficient smoothing part 51 in Fig. 5, 
30 wherein they are smoothed to obtain p smoothing coefficients s^k), s 2 (k), .... Sp.^k), s p (k). This smoothing operation is 

performed on the basis of Eq. (37). Of the smoothing coefficients, s^k), s 2 (k) Sp.^k) are supplied, as elements of 

the smoothing coefficient vector Sp^k), to the inner product calculation part 53 and an estimated transfer function cal- 
culation part 56. The smoothing coefficient s p (k) is fed to an estimated transfer function approximate value storage part 
55. 

35 [0035] The estimated transfer function approximate value storage part 55 updates the approximate value, using the 
smoothing coefficient s p (k) and the input signal vector x(k). That is, s p (k)x(k-p+1) is added to the approximate value 
z(k) stored until then and the added result is stored as z(k+1). This operation corresponds to the computation Eq. (40). 
[0036] Finally, in the estimated transfer calculation part 56 the smoothing coefficient vector s^k) and the input sig- 
nal vector x(k) are used to calculate Eq. (39) to obtain the estimated transfer function h (k+1). 

40 [0037] In the above-described operations, rough estimates of computational complexity involved in the respective 
parts are as follows: 



Correlation calculation part 52: < 


2p 


Inner product calculation part 53: 


P 


Convolution part 54: 


L 


Pre-filter coefficient vector calculation part 31: 


16p 


Pre-filter coefficient smoothing part 51: 


P 


Linear approximate value storage part 55: 


L 


Estimated transfer function calculation part 56: 


Lp 



55 

In the above, p-1 is regarded as nearly equal to p. The entire computational complexity is such as follows: 
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(L+20)p+2L 



(48) 



w 



15 



[0038] Now, attention should be paid to the following points. With the conventional transfer function estimation 
method depicted in Fig. 2. the estimated transfer function h (k+1) corresponding to the transfer function is calculated at 
every time and is used to synthesize the estimated value y (k) of the unknown system output In contrast thereto, 
according to the present embodiment which utilizes the approximate value z(k) of the estimated transfer function, the 
estimated value y (k) of the unknown system output can be synthesized without the need of calculating, directly, the 
estimated transfer function h (k+1) as shown in Fig. 5. Once the estimated value y (k) is obtained, the above-described 
operations can be performed. Accordingly, there is no need of calculating the estimated transfer function at every time 
in the cases where the estimated value h (k+1) of the transfer function, obtained as the result of the long-time estimat- 
ing operation, is needed and where the estimated value y (k) of the unknown system output is needed rather than the 
estimated result of the transfer function (for example, in the case of the estimation of characteristics of a time-invariant 
system or in an acoustic echo canceller). 

[0039] On this account, the overall computational complexity at every time, except the computational complexity in 
the estimated transfer function calculation part 56, is as follows: 



20p + 2L 



(49) 



20 



25 



30 



As referred to previously with respect to the prior art, the number of operations needed in the past is about 12000 when 
the tap number L of the filter is 500 and the projection order p is 20. From Eq. (49), however, the number of operations 
in the present invention is 1380; that is, the computational complexity is reduced to about 1/8 that in the prior art. 
[0040] As described above, the present invention permits substantial reduction of the computational complexity 
involved in the conventional projection scheme by the recursive synthesis of the pre-filter coefficient, storage of approx- 
imate values of the estimated transfer function vector and smoothing of the pre-filtering coefficients. 
[0041] While in the above the step size \i has been described to be handled as a scalar, it may also be provided as 
uA by use of a diagonal matrix A. For example, in the case where the energy of an impulse response of an unknown 
system decays exponentially, the transfer function estimation speed may sometimes be increased by arranging the ele- 
ments of the diagonal matrix A as a progression which decays exponentially, as shown below. 



A = diag(a, ay, ay 



,aY L " 1 ), (0<Y<D 



(50) 



In this instance, the input to the linear prediction part is multiplied by a ratio y 1/2 and R p (k) Eq. (14) is redefined by the 
following equation. 



35 



40 



R p (k) = 



x(k) A 
i <k-l) T 



x oc-p+ir 



A[x(k)x<k-1) ...*{**£=*) (51) 



45 This correction permits the application of the present invention described above. 

[0042] Next, a description will be given of examples of application of the transfer function estimating device of the 
present invention. 

[0043] A first example of application is the measurements of transfer functions of pieces of acoustic equipment. Fig. 
6 illustrates a system for measuring the transfer function from a loudspeaker to a microphone. In Fig. 6, reference 

so numeral 121 denotes a loudspeaker, 122 a microphone and 1 1 a transfer function estimating device. The output signal 
y(k) of the microphone 122 is a signal that has the characteristics of the loudspeaker 121 added to the input signal x(k). 
Regarding the loudspeaker 121 (including an acoustic path and the microphone as well) as an unknown system, the 
illustrated system is the same as that shown in Fig. 1, and by connecting the transfer function estimating device 11 of 
the present invention to the input of the loudspeaker 121 and the output of the microphone 122 as shown in Fig. 6, the 

55 transfer function h(k) of the loudspeaker can be estimated as the filter coefficient h (k) of the FIR filter. The measure- 
ment is usually performed in an anechoic chamber to avoid the influence of room acoustical characteristics. 
[0044] A second example is an acoustic echo canceller which prevents howling or echo in a loudspeaker commu- 
nication system such as a TV conference system or visual telephone. Fig. 7 is a diagram for explaining the acoustic 
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echo canceller. In Fig. 7, reference numeral 121 denotes a receiving loudspeaker, 122 a transmitting microphone, 123 
a room acoustic system and 20 an estimated echo generating part which is identical in construction with the estimated 
signal generating part 20 embodying the present invention in Fig. 2. The transfer function estimating device 1 1 of this 
embodiment operates as an acoustic echo canceller. In a hands-free communication system using the loudspeaker and 
5 the microphone, the other party's voice emanating from the receiving loudspeaker 121 is received by the transmitting 
loudspeaker via the room acoustic system of the transfer function h(k). The received signal y(k) is returned to the other 
party and reproduced there. At the other party side, the transmitted voice is sent back and reproduced; this phenome- 
non is called an acoustic echo and disturbs comfortable communication. 

[0045] The estimated echo generating part 20 of the acoustic echo canceller 11 estimates the room acoustical 
io characteristics h(k) including characteristics of the loudspeaker and generates an estimated echo signal y (k) based on 
the characteristics h (k) estimated as those of the input signal x(k). The subtracter 24 subtracts the estimated echo sig- 
nal y (k) from the received signal y(k). When the estimation is performed well, the echo canceller 1 1 operates so that it 
minimizes the power of the error signal e(k) and hence makes the estimated echo signal y (k) nearly equal to the 
received signal y(k), substantially reducing the acoustic echo. 
15 [0046] Comparison of the Fig. 7 system with the Fig. 2 system reveals that the transfer function estimating device 
of the present invention can be used directly as an echo canceller. The room acoustic system 123 in Fig. 7 corresponds 
to the unknown system 12 in Fig. 2 and the estimated echo y (k) and the transmitted signal e(k) in Fig. 7 correspond to 
the output y(k) from the convolution part 23 and the error signal e(k) in Fig. 2, respectively. 

[0047] A third example is noise control. Fig. 8 illustrates the principles of noise control. In Fig. 8, reference numeral 
20 131 denotes a noise source, 132 a noise transfer path expressed by the transfer function h(k), 133 a microphone for 
observation, 134 a noise monitor microphone, 20 an estimated noise generating part, 136 a phase inverter and 137 a 
loudspeaker. The purpose of noise control is to cancel noise that is observed via the noise transfer path of the transfer 
characteristics h(k) from the noise source 131, by generating negative estimated noise (letting a sound pressure be rep- 
resented by y(k), a sound pressure represented by -y(k) is called a negative sound with respect to y(k)) from the loud- 
25 speaker 137. 

[0048] To attain this object, the transfer characteristics h(k) of the noise transfer path 132 is estimated in the esti- 
mated noise generating part 20 of the same construction as that of the estimated signal generating part 20 in Fig. 2. 
That is, noise is detected by the monitor microphone 134 in the vicinity of the noise source 131 and provided as the 
input signal x(k) to the estimated noise generating part 20, which generates an estimated value y (k) of the noise signal 

30 y(k) at the observation point (i.e. the microphone 133). The estimated noise y (k) is polarity inverted by the phase 
inverter 136 into a signal -y(k). Assuming, for the sake of brevity, that the loudspeaker characteristics are negligible, the 
signal -y (k) from the loudspeaker 137 is a combination of the noise signal y(k) and the sound pressure in the micro- 
phone 133 at the observation point and the error signal e(k) is provided at the output of the microphone 133. In this 
instance, when the noise transfer path is estimated well, the signal y (k) becomes similar to the noise y(k) from the 

35 noise source 131 and the noise y(k) is cancelled by the combined sound pressure -y (k) from the loudspeaker 137. 
[0049] The microphone 133 in Fig. 8 corresponds to the subtracter 24 in Fig. 2; accordingly, the systems of Figs. 8 
and 2 differ only in whether the error signal generated outside the estimating device 1 1 is provided thereto or the error 
signal is calculated in the estimating device 1 1 . The principles of the present invention are applicable to the noise con- 
trol device of Fig. 8. 

40 [0050] In the measurement of the transfer characteristics of the loudspeaker shown in Fig. 6, the loudspeaker char- 
acteristics do not change with time; hence, there is no need of learning measured results (transfer characteristics) in 
the course of measurement and the transfer characteristics need only to be made known as the result of measurement 
conducted after a certain period of time. In the acoustic echo canceller shown in Fig. 7, the room transfer function varies 
with time in response to movement of audiences or opening and closing of a door, for instance. As is evident from Fig. 

45 7, however, the purpose of the acoustic echo canceller is attained by obtaining the estimated value y (k) of the output 
y(k) of the room transfer function (the unknown system). Therefore, the estimated value h (k) of the room transfer func- 
tion itself is unnecessary in the acoustic echo canceller. Also in the noise control device of Fig. 8, its purpose can be 
attained by obtaining the estimated value y (k) of the unknown system and no estimated value of the room transfer func- 
tion itself is needed. 

50 [0051] As will be appreciated from the above, the transfer function estimating device of the present invention, when 
applied as described above, permits substantial reduction of the computational complexity as compared with the prior 
art. Furthermore, the present invention has a basic feature of minimizing the power of the error signal e(k) in such a 
system configuration as shown in Fig. 2. Thus, the present invention is applicable to all instances in which problems to 
be solved can be modelled as the error signal power minimizing problem shown in Fig. 2. 

55 [0052] In the embodiments described above, it will be effective to select the elements of the diagonal matrix A of 
Eq. (50) so that they provide the same attenuation as the room reverberation. 

[0053] Finally, a description will be given of experimental results of the acoustic echo canceller according to the 
present invention. Fig. 9 shows what are called learning curves, the ordinate representing echo return loss enhance- 
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ment (hereinafter referred to as ERLE) and the abscissa representing time. As the estimation of the room transfer func- 
tion proceeds with the lapse of time, the ERLE increases. In Fig. 9, the curve 21 1 is a learning curve in the case of the 
projection order p being 2, the curve 212 a learning curve in the case of the projection order p being 8 and the curve 
213 a learning curve in the case of the projection order p being 32. It will be understood, from Fig. 9, that the larger the 

5 projection order, the higher the convergence speed (the ERLE increases in a short period of time). 

[0054] Fig. 10 is a graph showing the relationship between the projection order p and the required computational 
complexity, the abscissa representing the projection order p and the ordinate the number of multiply-add operations 
(including add operations as well). In Fig. 10, the curve 221 indicates the computational complexity in the prior art and 
the curve 222 the computational complexity when the present invention was used. The tap number L of the FIR filter for 

10 use in the estimation of the transfer function was 500. From Fig. 10 it will be seen that the computational complexity 
involved in the present invention is much smaller than in the case of the prior art when the projection order p is large. 
[0055] As described above, the present invention permits substantial reduction of the number of operations neces- 
sary for the estimation of the transfer function or output of an unknown system by the projection scheme. In concrete 
terms, letting the tap length of the FIR filter which represents the unknown system and the projection order be repre- 

15 sented by L and p, respectively, the prior art requires multiply-add operations p 3 /6+(p+1)L , that is, the computational 
complexity is proportional to p 3 , whereas according to the present invention the computation complexity can be reduced 
down to 20p+2L. 

[0056] Such a reduction of computational complexity allows corresponding reduction in the scale of hardware, and 
hence significantly contributes to the downsizing of the device and reduction of its cost. Moreover, when the hardware 
20 is held on the same scale, the tap length L and the projection order can be chosen larger than in the past - this speeds 
up the estimating operation and increases the estimation accuracy. Besides, when the present invention is embodied 
by a computer, the operation time can be greatly shortened. 

Claims 

25 

1 . A method of adaptively estimating the transfer function of an unknown technical system by use of a projection algo- 
rithm, said method comprising: 

(i) applying an input physical signal x(k) to said unknown system (12), k representing a discrete time; 
30 (ii) detecting an output physical signal y(k) from said unknown system (12); 

(iii) estimating the transfer function using said input signal x(k) and said output signal y(k) and generating an 
output signal y (k) of an estimated system having said estimated transfer function h (k), wherein said estimated 
system is implemented as a digital filter of tap number L; 

(iv) obtaining, as an error signal e(k)=y(k)-y(k) , the difference between the output signal y (k) of said esti- 
35 mated system and the output signal y(k) of said unknown system (1 2); 

(v) calculating a correcting vector 8h (k) for the estimated transfer function from input signal vectors x(k), x(k- 

1), .... x(k-p + 1) and error signals e(k), e(k-1) e(k-p + 1) , where p is an integer equal to or greater than 2 

and x(k) is defined as x(k)=[x(k), x(k-1) x(k-L+ 1)] ; and 

(vi) using said correcting vector 5h (k) and a predetermined correcting step size u to repeatedly correct said 
40 estimated transfer function h (k) at each time k by the following equation: 

h(k+1) = h(k) + u8h(k) 

so that said error signal e(k) approaches zero; 
45 characterized in that step (v) comprises the steps of: 

(v-1) calculating a pre-filter coefficient vector g(k) by solving the following simultaneous linear equation 
with p unknowns: 

so R p (k)g(k)=e(k), 

wherein e(k) represents the vector of said error signal e(k) and is defined as 
e(k)=[e(k), (1-u)e(k-1) (l-uj^efk-p+l)] 1 . and 

55 

R p (k) is a covariance matrix of said input signal x(k) defined as 

R p (k) = Wk). x(k-1) x(k-p + 1)] T [x(k), x(k-1),..., x(k-p + 1)1; and 



15 



EP 0 637 803 B1 



(v-2) using said pre-filter coefficient vector 9(k) to calculate said correcting vector 6h (k) by the following 



wherein x(k) represents the vector of said input signal x(k); 
wherein step (v-1) comprises the steps of: 

(v-1-a) calculating a forward linear prediction coefficient vector a(k) of said input signal x(k), the sum of its 
a posteriori prediction-error squares F(k), a backward linear prediction coefficient vector b(k) of said input 
signal x(k) and the sum of its a posteriori prediction-error squares B(k), and 

(v-1-b) letting a pre-filter deriving coefficient vector be represented by f(k), obtaining said pre-filter coeffi- 
cient vector g(k) by a recursion formula composed of the following first and second equations: 



A method of adaptively estimating the transfer function of an unknown technical system by use of a projection algo- 
rithm, said method comprising: 

(i) applying an input physical signal x(k) to said unknown system (12), k representing a discrete time; 

(ii) detecting an output physical signal y(k) from said unknown system (12); 

(iii) estimating the transfer function using said input signal x(k) and said output signal y(k) and generating an 
output signal y (k) of an estimated system having said estimated transfer function h (k), wherein said estimated 
system is implemented as a digital filter of tap number L; 

(iv) obtaining, as an error signal e(k)=y(k)-y(k) , the difference between the output signal y (k) of said esti- 
mated system and the output signal y(k) of said unknown system (12); 

(v) repeating the calculation of the estimated transfer function at each time k so that said error signal e(k) 
approaches zero; 

characterized in that 
step (v) comprises steps of: 

(v-1) calculating a pre-filter coefficient vector g(k) by solving the following simultaneous linear equation 
with p unknowns: 



equation: 



8h(k) =[x(k),x(k-1) x(k-p+1)]g(k), 





R p (k)g(k)=e(k), 



wherein e(k) represents the vector of said error signal e(k) and is defined as . 



e(k)=[e(k), (1-u)e(k-1) (l-nJ^efk-p+l)] 1 , and 



R p (k) is a covariance matrix of said input signal x(k) defined as 



R p (k) - Mk), x(k-1) x(k-p+1)] T (x(k), x(k-1) x(k-p+1)]; 



16 



EP 0 637 803 B1 



(v-2) smoothing pre-fitter coefficients gj(k) which are elements of said pre-filter coefficient vector g(k) by 
the following equations: 

s i (k)=s M (k-1) + ug j (k) for2^p 
^g^k) for i=1 



to obtain a smoothing coefficient Sj(k); 

(v-3) obtaining an approximate estimated transfer function z(k + 1) as an approximation of said estimated 
transfer function using the smoothing coefficient s p (k) based on the following equation 

z(k + 1)=z(k) + s p (k)x(k-p+1); 

(v-4) calculating the convolution x(k) T z(k) of said approximate estimated transfer function z(k) with an input 
signal vector x(k) defined as x(k)=[x(k), x(k-1),... f x(k-L + 1)] ; 

(v-5) calculating the inner product s ^ (k-1) r p _ 1 (k) , setting the vector of said smoothing coefficient srfk) 
and a correlation vector r^k) of said input signal to 



and 



Sp.^k-IJ^s^k-IJ.s^M) s^k-l)] 1 . 

r p . 1 (k)=[x(k) T x(k-1), x(k) T x(k-2) x(k) T x(k-p+1)] T t 



respectively; and 

(v-6) obtaining the sum of said convolution result x(k) T z(k) and said inner product as said estimated output 
.signal y (k). 

The method of claim 2, wherein step (v-1) comprises: 

(a) calculating a forward linear prediction coefficient vector a(k) of said input signal x(k), the sum of its a pos- 
teriori prediction-error squares F(k), a backward linear prediction coefficient vector b(k) of said input signal x(k) 
and the sum of its a posteriori prediction-error squares B(k), and 

(b) letting a pre-filter deriving coefficient vector be represented by f(k), obtaining said pre-filter coefficient vector 
g(k) by a recursion formula composed of the following first and second equations: 



»-* i, "j, T *-v -< 



The method of claim 1 or 3, wherein, letting the last element of said pre-filter coefficient vector g(k-1) be repre- 
sented by g p (k-1), said pre-filter coefficient vector g(k) is calculated using the following equation which is a modified 
version of said second equation: 




= g(k-l)-g p (k-i)b(k-1) 



The method of claim 1 or 3, wherein, letting predetermined two non-negative numbers be represented by 6 F (k) and 
8 B (k-1), said pre-filter coefficient vector g(k) is calculated by the following equations which are modified versions of 
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said first and second equations: 



r 0 1 a(kTe(k) 
|f(k-1)J F(k)+8 F (k) 



[V]-*-'- 



b(k-1) T e(k-1) 
B(k-1)+8 B (k-1) 



T b(k-1) 



The method of claim 2 or 3 wherein step (v-5) comprises calculating said correlation vector r^k) of said input sig- 



A device for adaptive estimation of the transfer function of an unknown technical system by use of a projection algo- 
rithm, said device comprising: 

means for applying an input physical signal x(k) to said unknown system (121 , k representing a discrete time; 
means for detecting an output physical signal y(k) from said unknown system (12); 

means (1 1) for estimating the transfer function using said input signal x(k) and said output signal y(k) and gen- 
erating an output signal y (k) of an estimated system (23) having said estimated transfer function h (k), wherein 
said estimated system is implemented as a digital filter of tap number L; 

means (24) for obtaining, as an error signal e(k)=y(k)-y(k) , the difference between the output signal y (k) of 
said estimated system and the output signal y(k) of said unknown system (12); 

means (21) for calculating a correcting vector 8h (k) for the estimated transfer function from input signal vectors 

x(k), x(k-1), ... , x(k-p+1) and error signals e(k), e(k-1) e(k-p+1 ), where p is an integer equal to or greater 

than 2 and x(k) is defined as x(k)=[x(k), x(k-1) x(k- L+1)] ; and 

means (22) for repeatedly correcting, at each time k using said correcting vector 5h (k) and a predetermined 
correcting step size u, said estimated transfer function h (k) by the following equation: 



so that said error signal e(k) approaches zero; 

characterized in that said means (21) for calculating the correcting vector 6h (k) comprises: 

- means (31) for calculating a pre-filter coefficient vector g(k) by solving the following simultaneous linear 
equation with p unknowns: 



nal x(k) by 



'p-i(k) = r p -i(k-1)-x(k-L)x p . 1 (k-L)+x(k)x p . l (k) 
where x^k) represents the vector of said input signal x(k) expressed as: 

Xp^WWk-1). x(k-2) x(k-p+1)] T . 



h(k+1) = h(k) + u5h(k) 



R p (k)g(k)=e(k), 

wherein e(k) represents the vector of said error signal elk) and is defined as 



e(k)=[e(k), (1-u)e(k-1) (1^) p * 1 e(k-p+1)] T , and 

R p (k) is a covariance matrix of said input signal x(k) defined as 

R p (k)=[x(k), x(k-1) x(k-p+1)] T [x(k), x(k-1) x(k-p+1)]; and 



means (32) for calculating, using said pre-filter coefficient vector g(k), said correcting vector 6h (k) by the 
following equation: 
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6h(k) =[x(k),x(k-1) x(k-p+1)]g(k), 



where x(k) represents the vector of said input signal x(k); and 

said means (31) for calculating said pre-filter coefficient vector g(k) comprises: 

means (41) for calculating a forward linear prediction coefficient vector a(k) of said input signal x(k), the 
sum of its a posteriori prediction-error squares F(k), a backward linear prediction coefficient vector b(k) of 
said input signal x(k) and the sum of its a posteriori prediction-error squares B(k), and 
means (42, 43. 44) for obtaining said pre-filter coefficient vector g(k) by a recursion formula composed of 
the following first and second equations: 



where f(k) represents a pre-filter deriving coefficient vector. 

A device for adaptive estimation of the transfer function of an unknown technical system by use of a projection algo- 
rithm, said device comprising: 

means for applying an input physical signal x(k) to said unknown system (12), k representing a 
discrete time; 

means for detecting an output technical signal g(k) from said unknown system (12); 

means (1 1) for estimating the transfer function using said input signal x(k) and said output signal y(k) and gen- 
erating an output signal y (k) of an estimated system (23) having said estimated transfer function h (k), wherein 
said estimated system is implemented as a digital filter of tap number L; 

means (24) for obtaining, as an error signal e(k) = y(k)-y(k) , the difference between the output signal y (k) of 
said estimated system and the output signal y(k) of said unknown system (12); and 

means (20) for calculating the estimated transfer function at each time k so that said error signal e(k) 
approaches zero; 

characterized In that said means (20) for calculating the estimated transfer function comprises: 

means (31) for calculating a pre-filter coefficient vector g(k) by solving the following simultaneous linear 
equation with p unknowns: 



means (51 ) for smoothing pre-filter coefficients g^k), which are elements of said pre-filter coefficient vector 
g(k), by the following equations: 




R p (k)g(k) = e(k), 



wherein e(k) represents the vector of said error signal e(k) and is defined as 



e(k)=[e(k), (1-u)e(k-1) (l-uj^efk-p+l)] 1 , and 



R p (k) is a covariance matrix of said input signal x(k) defined as 



R p (k) = [x(k), x(k-1) x(k-p + 1)] T (x(k), x(k-1) x(k-p+1)]; 



Sj (k)=s M (k-1)+u gi (k) 
= ug 1 (k) for i=1 



for2i£i2Sp 



to obtain a smoothing coefficient Sj(k); and 
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- means (55) for obtaining an approximate estimated transfer function z(k + 1) as an approximation of said 
estimated transfer function using the smoothing coefficient s p (k) based on the following equation 

z(k ♦ 1) = z(k) + s p (k)x(k-p + 1); and 

5 

said estimated system (23) comprises: 

- means (54) for calculating the convolution x(k) T z(k) of said approximate estimated transfer function z(k) 
with an input signal vector x(k) defined as x(k)=[x(k), x(k-1),..., x(k- L+1)] T ; 

10 - means (53) for calculating the inner product s ^(k-1) r p . 1 (k) , setting the vector of said smoothing coef- 

ficient Sj(k) and a correlation vector r^k) of said input signal to 

s p.1 (k-1 )=[s , (k-1), s 2 (k-1 ) s M (k-1 )] T 

15 and 

r p . 1 (k)=[x(k) T x(k-1). x(k) T x(k-2) x(k) T x(k-P + 1)] T 

respectively; and 

20 - means (57) for obtaining the sum of said convolution result x(k) T z(k) and said inner product as said esti- 

mated output signal y (k). 

9. The device of claim 8, wherein said means (31 ) for solving said simultaneous linear equation with p unknowns com- 
prises: 

25 

means (41) for calculating a forward linear prediction coefficient vector a(k) of said input signal x(k), the sum of 
its a posteriori prediction-error squares F(k), a backward linear prediction coefficient vector b(k) of said input 
signal x(k) and the sum of its a posteriori prediction-error squares B(k), and 

means (42, 43, 44) for obtaining said pre-filter coefficient vector g(k) by a recursion formula composed of the 
30 following first and second equations: 



35 



40 



where f(k) represents a pre-filter deriving coefficient vector. 

10. The device of claim 8 or 9, wherein said inner product calculating means includes correlation calculating means for 
calculating said correlation vector ^(k) of said input signal x(k) by 

Vl< k > = r p . 1 (k-1)-x(k-L)x p . 1 (k-L) + x(k)x p . 1 (k), 

45 

where x p _ 1 (k) represents the vector of said input signal x(k) expressed as: 

x p . 1 (k) = [x(k-1), x(k-2) x(k-p+1)] T . 

so 11. The device of claim 7 or 9, wherein said pre-filter deriving coefficient vector correcting means is means which, let- 
ting the last element of said pre-filter coefficient vector g(k-1) be represented by g p (k-1), calculates said pre-filter 
coefficient vector g(k) by the following equation which is a modified version of said second equation: 



55 



[ ,0t 0 ' 1) ]=g(k-D-gp(k-i)b{k-i) 
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12. The device of claim 7 or 9, wherein, letting predetermined non-negative numbers be represented by 5p(k) and ^(k- 
1), respectively, said pre-filter coefficient vector correcting means and said pre-filter deriving coefficient vector cor- 
recting means calculate said pre-filter coefficient vector g(k) by the following equations, respectively, which are 
modified versions of said first and second equations: 



13. The device of claim 7 or 10, further comprising subtracting means for providing, as said error signal e(k), the differ- 
ence y(k)-y (k) between said output signal y(k) of said unknown system (12) and said estimated signal y (k) which 
is the output of said estimated system. 

Patentanspruche 

1. Veriahren zum adaptiven Schatzen der Ubertragungsfunktion eines unbekannten technischen Systems durch Ver- 
wendung eines Projektionsalgorithmus, wobei das Verfahren umfafit: 

(i) Anlegen eines physikalischen Eingangssignals x(k) an das unbekannte System (12), wobei k eine diskrete 
Zeit reprasentiert; 

(it) Erfassen eines physikalischen Ausgangssignals y(k) aus dem unbekannten System (12); 
(Hi) Schatzen der Ubertragungsfunktion unter Verwendung des Eingangssignals x(k) und des Ausgangssignals 
y(k) und Erzeugen eines Ausgangssignals y (k) eines Schatzsystems mit der geschatzten Ubertragungsfunk- 
tion h (k), wobei das Schatzsystem als ein digitales Filter mit einer Anzahl L an Anzapfungen implementiert ist; 

(iv) Ermitteln, als ein Fehlersignal e(k) = y(k)-y(k) , der Differenz zwischen dem Ausgangssignal y (k) des 
Schatzsystems und dem Ausgangssignal y(k) des unbekannten Systems (12); 

(v) Berechnen eines Korrigiervektors 5h (k) fur die geschatzte Ubertragungsfunktion aus Eingangssignalvek- 

toren x(k), x(k-1) x(k-p + 1) und Fehlersignalen e(k), e(k-1), ... t e(k-p + 1) , wobei p eine ganze Zahl grdBer 

oder gleich 2 ist und x(k) definiert ist als x(k) = [x(k), x(k-1) x(k-L + 1)] T ; und 

(vi) Verwenden des Korrigiervektors 5h (k) und einer vorbestimmten Korrigierschrittgrofie \i, urn wiederholt die 
geschatzte Ubertragungsfunktion h (k) zu jedem Zeitpunkt k durch die fblgende Gleichung zu korrigieren: 



so daft sich das Fehlersignal e(k) Null annahert; 

dadurch gekennzeichnet, dad Schritt (v) fblgende Schritte umfafit: 

(v-1) Berechnen eines Vorfilterkoeffizientenvektors g(k) durch Losen der fblgenden simultanen linearen 
Gleichung mit p Unbekannten: 




h(k+1)=h(k)+u5h(k) 



R p (k)g(k)=e(k), 



wobei e(k) den Vektor des Fehlersignals e(k) reprasentiert und definiert ist als 



e(k)=[e(k), (1-u)e(k-1) (l-uj^efk-p+l)] 1 , und 



R p (k) eine Kovarianzmatrix des Eingangssignals x(k) ist und definiert ist als 



R p (k)=[x(k), x(k-1) x(k-p + 1)] T [x(k), x(k-1) x(k-p + 1)|; und 



(v-2) Verwenden des Vorfilterkoeffizientenvektors g(k) zum Berechnen des Korrigiervektors 8h (k) durch 
die folgende Gleichung: 
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5h(kMx(k), X (k-1) x(k-p+1)]g(k), 



wobei x(k) den Vektor des Eingangssignals x(k) reprasentiert; 
wobei Schritt (v-1) folgende Schritte umfa&t: 

(v-1-a) Berechnen eines Vorvvartslinearvorhersagekoeffizientenvektors a(k) des Eingangssignals x(k), der 
Summe seiner a-posteriori-Vorhersagefehlerquadrate F(k), eines Ruckwartslinearvorhersagekoef- 
fizientenvektors b(k) des Eingangssignals x(k) und der Summe seiner a-posteriori-Vorhersagefehlerqua- 
drate B(k), und 

(v-1-b) wenn ein Vorfirterableitungskoeffizientenvektor durch f(k) reprasentiert sei, Ermitteln des Vorfifter- 
koeffizientenvektors g(k) durch eine Rekursionsformel, die sich aus der folgenden ersten und zweiten Glei- 
chung zusammensetzt: 



Verfahren zum adaptiven Schatzen der Obertragungsfunktion eines unbekannten technischen Systems durch Ver- 
wendung eines Projektionsalgorithmus, wobei das Verfahren umfaBt: 

(i) Anlegen eines physikalischen Eingangssignals x(k) an das unbekannte System (12), wobei k eine diskrete 
Zeit reprasentiert; 

(ii) Erfassen eines physikalischen Ausgangssignals y(k) aus dem unbekannten System (12); 

(iii) Schatzen der Obertragungsfunktion unter Verwendung des Eingangssignals x(k) und des Ausgangssignals 
y(k) und Erzeugen eines Ausgangssignals y (k) eines Schatzsystems mit der geschatzten Ubertragungsfunk- 
tion h (k), wobei das Schatzsystem als ein digitales Filter mit einer Anzahl L an Anzapfungen implementiert ist; 

(iv) Ermitteln, als ein Fehlersignal e(k)=y(k)-y(k) , der Differenz zwischen dem Ausgangssignal y (k) des 
Schatzsystems und dem Ausgangssignal y(k) des unbekannten Systems (12); 

(v) Wiederholen der Berechnung der geschatzten Obertragungsfunktion zu jedem Zeitpunkt k so, dad sich das 
Fehlersignal e(k) Null annahert; 

dadurch gekennzeichnet, dad 
Schritt (v) folgende Schritte umfafM: 

(v-1) Berechnen eines Vorfilterkoeffizientenvektors g(k) durch Losen der folgenden simultanen linearen 
Gleichung mit p Unbekannten: 



R p (k) eine Kovarianzmatrix des Eingangssignals x(k) ist und definiert ist als 

R p (k) = [x(k), x(k-1) x(k-p + 1)] T [x(k), x(k-1) x(k-p + 1)J; 

(v-2) Glatten von Vorfilterkoeffizienten g(k), die Elemente des Vorfilterkoeffizientenvektors g(k) sind, durch 
die folgenden Gleichungen: 





R p (k)g(k)=e(k), 



wobei e(k) den Vektor des Fehlersignals e(k) reprasentiert und definiert ist als 



e(k) =[e(k), (1-n)e(k-1) (1-u) p " 1 e(k-p + 1)] T , und 



8j(k)=3 M (k-1)+u g| (k) 
=u. gi (k) furi=1 



fur2<i^p 
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um einen Glattungskoeffizienten Sj(k) zu ermitteln; 

(v«3) Ermitteln einer ungefahren geschatzten Obertragungsfunktion z(k + 1) als einer Naherung der 
geschatzten Obertragungsfunktion unter Verwendung des Glattungskoeffizienten s p (k) auf der Basis der 
folgenden Gleichung 

z(k + 1)=z(k) + s p (k)x(k-p+1); 

(v-4) Berechnen der Faltung x(k) T z(k) der ungefahren geschatzten Obertragungsfunktion z(k) mit einem 
Eingangssignatvektor x(k), der definiert ist als x(k) = [x(k), x(k-1), .... x(k-L + 1 ) ; 
(v-5) Berechnen des inneren Produkts Sp.^k-1) r^k), wobei der Vektor des Glattungskoeffizienten 
Sj(k) und ein Korrelationsvektor r^k) des Eingangssignals auf 

s p . 1 (k-1)=[s 1 (k-1),s 2 (k.1) s p . l (k.1)] T , 

bzw. 

r M (k)=[x(k) T x(k-1), x(k) T x(k-2) x(k) T x(k-p + 1)] T t 

gesetzt werden; und 

(v-6) Ermitteln der Summe des Faltungsergebnisses x(k) T z(k) und des inneren Produkts als das 
geschatzte Ausgangssignal y (k). 

Verfahren nach Anspruch 2, bei dem Schritt (v-1) umfa&t: 

(a) Berechnen eines VorwartsHnearvorhersagekoeffizientenvektors a(k) des Eingangssignals x(k), der Summe 
seiner a-posteriori-Vorhersagefehlerquadrate F(k), eines Ruckwartslinearvorhersagekoeffizientenvektors b(k) 
des Eingangssignals x(k) und der Summe seiner a-posteriori-Vomersagefehlerquadrate B(k), und 

(b) wenn ein Vorfilterableitungskoeffizientenvektor durch f(k) reprasentiert sei, Ermitteln des Vorfilterkoeffizlen- 
tenvektors g(k) durch eine Rekursionsformel, die sich aus der folgenden ersten und zweiten Gleichung zusam- 
mensetzt: 



[V] -•*■>• 



Verfahren nach Anspruch 1 Oder 3, bei dem, wenn das letzte Element des Vorfilterkoeffizientenvektors g(k-1) als 
durch g p (k-1) reprasentiert angenommen wird, der Vorfilterkoeffizientenvektor g(k) unter Verwendung der folgen- 
den Gleichung berechnet wird, die eine modifizierte Version der zweiten Gleichung ist: 




= gOc-1) -g p <k-u b(k-D 



Verfahren nach Anspruch 1 Oder 3, bei dem, wenn zwei vorbestimmte nicht-negative Zahlen durch 5p(k) und 5e(k- 
1) reprasentiert seien, der Vorfilterkoeffizientenvektor g(k) durch die folgenden Gleichungen berechnet wird, die 
modifizierte Versionen der ersten und der zweiten Gleichung sind: 
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[ o J 91 B(k-l)+5B(k-D 



15 

6. Verfahren nach Anspruch 2 Oder 3, bei dem Schrift (v-5) das Berechnen des Korrelationsvektors r^k) des Ein- 
gangssignals x(k) durch 



20 



f p-iM = r p . 1 (k-1)-x(k-L)x p . 1 (k-L) + x(k)x p . 1 (k) 

umfaRt, wobei x^^k) den Vektor des Eingangssignals x(k) reprasentiert und ausgedruckt ist als: 

x M (k)=[x(k-1) f x(k-2) x(k-p+1)] T . 

2$ 7. Vorrichtung zum adaptiven Schatzen der Ubertragungsfunktion eines unbekannten technischen Systems durch 
Verwendung eines Projektionsalgorithmus, wobei die Vorrichtung umfaBt: 

Mittel zum Anlegen eines physikalischen Eingangssignals x(k) an das unbekannte System (12), wobei k eine 
diskrete Zeit reprasentiert; 

30 Mittel zum Erfassen eines physikalischen Ausgangssignals y(k) aus dem unbekannten System (12); 

Mittel (11) zum Schatzen der Ubertragungsfunktion unter Verwendung des Eingangssignals x(k) und des 
Ausgangssignals y(k) und Erzeugen eines Ausgangssignals y (k) eines Schatzsystems (23) mit der geschatz- 
ten Ubertragungsfunktion h (k), wobei das Schatzsystem als ein digitales Filter mit einer Anzahl L an Anzap- 
fungen implementiert ist; 

35 Mittel (24) zum Ermitteln, als ein Fehlersignal e(k)=y(k)-y(k) , der Differenz zwischen dem Ausgangssignal 

y (k) des Schatzsystems und dem Ausgangssignal y(k) des unbekannten Systems (12); 
Mittel (21) zum Berechnen eines Korrigiervektors 8h(k) fur die geschatzte Ubertragungsfunktion aus Ein- 

gangssignalvektoren x(k), x(k-1) x(k-p + 1) und Fehlersignalen e(k), e(k-1), .... e(k-p + 1) , wobei p eine 

ganze Zahl grdlier oder gleich 2 ist und x(k) definiert 1st als x(k) = [x(k) f x(k-1) x(k-L + 1)] ; und 

40 Mittel (22) zum wiederholten Korrigieren, unter Verwendung des Korrigiervektors 8h (k) und einer vorbestimm- 

ten Korrigierschringrolie u, der geschatzten Obertragungsfunktion h (k) zu jedem Zeitpunkt k durch die fol- 
gende Gleichung: 



45 



h(k+ 1)=h(k)+u5h(k) 
so daft sich das Fehlersignal e(k) Null annahert; 

dadurch gekennzelchnet dali die Mittel (21) zum Berechnen des Korrigiervektors 8h (k) umfassen: 

Mittel (31) zum Berechnen eines Vorfilterkoeffizientenvektors g(k) durch Losen der folgenden simultanen 
so linearen Gleichung mit p Unbekannten: 

R p (k)g(k)=e(k), 



55 



wobei e(k) den Vektor des Fehlersignals e(k) reprasentiert und definiert ist als 

e(k)=[e(k), (1-u)e(k-1), .... (l-^^etk-p+l)] 1 , und 
Rp(k) eine Kovarianzmatrix des Eingangssignals x(k) ist und definiert ist als 
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R p (k)=Mk). x(k-1) x(k-p + 1)] T [x(k), x(k-1) x(k-p + 1)]; und 

Mittel (32) zum Berechnen des Korrigiervektors 8h (k) unter Verwendung des Vorfilterkoeffizientenvektors 
g(k) durch die fbigende Gleichung: 

8h(k)=[x(k), x(k-1) x(k-p+1)]g(k), 

wobei x(k) den Vektor des Eingangssignals x(k) reprasentiert; 

wobei die Mittel (31) zum Berechnen des Vorfilterkoeffizientenvektors g(k) umfaftt: 

Mittel (41) zum Berechnen eines Vorwartslinearvorhersagekoeffizientenvektors a(k) des Eingangssignals 

x(k) t der Summe seiner a-posteriori-Vorhersagefehlerquadrate F(k), eines Ruckwartslinearvorhersageko- 

effizientenvektors b(k) des Eingangssignals x(k) und der Summe seiner a-posteriori-Vorhersagefehlerqua- 

drate B(k), und 

Mittel (42. 43, 44) zum Ermitteln des Vorfilterkoeffizientenvektors g(k) durch eine Rekursionsformel, die 
sich aus der folgenden ersten und zweiten Gleichung zusammensetzt: 



['V']- 



d(K- i; 



wobei f(k) einen Vorfilterableitungskoeffizientenvektor reprasentiert. 

Vorrichtung zum adaptiven Schatzen der Ubertragungsfunktion eines unbekannten technischen Systems durch 
Verwendung eines Projektionsalgorithmus, wobei die Vorrichtung umfaftt: 

Mittel zum Anlegen eines physikalischen Eingangssignals x(k) an das unbekannte System (12), wobei k eine 
diskrete Zeit reprasentiert; 

Mittel zum Erfassen eines physikalischen Ausgangssignals y(k) aus dem unbekannten System (12); 
Mittel (11) zum Schatzen der Obertragungsfunktion unter Verwendung des Eingangssignals x(k) und des 
Ausgangssignals y(k) und Erzeugen eines Ausgangssignals y (k) eines Schatzsystems (23) mit der geschatz- 
ten Ubertragungsfunktion h (k), wobei das Schatzsystem als ein digitales Filter mit einer Anzahl L an Anzap- 
fungen implementiert ist; 

Mittel (24) zum Ermitteln, als ein Fehlersignal e(k) = y(k)- y(k) , der Differenz zwischen dem Ausgangssignal 
y (k) des Schatzsystems und dem Ausgangssignal y(k) des unbekannten Systems (12); und 
Mittel (20) zum Berechnen der geschatzten Ubertragungsfunktion zu jedem Zeitpunkt k so, daft sich das Feh- 
lersignal e(k) Null annahert; 

dadurch gekennzeichnet, daft Mittel (20) zum Berechnen der geschatzten Ubertragungsfunktion umfassen: 
Mittel (31) zum Berechnen eines Vorfilterkoeffizientenvektors g(k) durch Losen der folgenden simultanen linea- 
ren Gleichung mit p Unbekannten: 

R p (k)g(k)=e(k), 

wobei e(k) den Vektor des Fehlersignals e(k) reprasentiert und definiert ist als 

e(k) = [e(k), (1-^i)e(k-1), .... (l-jO^efk-p + 1)J T , und 

R p (k) eine Kovarianzmatrix des Eingangssignals x(k) ist und definiert ist als 

R p (k) = (x(k), x(k-1) x(k-p + 1)] T [x(k), x(k-1) x(k-p + 1)]; 

Mittel (51) zum Glatten von Vorfilterkoeffizienten gj(k), die Elemente des Vorfilterkoeffizientenvektors g(k) 
sind, durch die folgenden Gleichungen: 
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srfk^s^k-IJ + ugrfk) fur2*i*p 
=Mi(k) furi =1 



zum Ermitteln eines Glattungskoeffizienten Sj(k); und 

- Mittel (55) zum Ermitteln einer ungefahren geschitzten Obertragungsfunktion z(k + 1) als eine Naherung 
der geschatzten Obertragungsfunktion unter Verwendung des Glattungskoeffizienten s p (k) auf der Basis 
der fotgenden Gleichung 

z(k+1)=z(k) + s p (k)x(k-p+1); 

wobei das Schatzsystem (23) umfaBt: 

- Mittel (54) zum Berechnen der Faltung x(k) T z(k) der ungefahren geschatzten Obertragungsfunktion z(k) 
mit einem Eingangssignalvektor x(k), der definiert ist als x(k) = [x(k), x(k-1) x(k-L +1)] ; 

- Mittel (53) zum Berechnen des inneren Produkts s M (k-1) r ^ (k) , wobei der Vektor des Glattungskoef- 
fizienten Sj(k) und ein Korrelationsvektor r^(k) des Eingangssignals auf 

Sp^k-IHs^k-l), s 2 (k-1) s p . 1 (k.1)] T , 

bzw. 

r M (k) = [x(k) T x(k-1), x(k) T x(k-2) x(k) T x(k-P ♦ 1)] T , 

gesetzt werden; und 

- Mittel (57) zum Ermitteln der Summe des Faltungsergebnisses x(k) T z(k) und des inneren Produkts als das 
geschatzte Ausgangssignal y (k). 

9. Vorrichtung nach Anspruch 8, bei der die Mittel (31 ) zum Losen der simultanen linearen Gleichung mit p Unbekann- 
ten umfaftt: 

- Mittel (41) zum Berechnen eines Vorwartslinearvorhersagekoeffizientenvektors a(k) des Eingangssignals x(k), 
der Summe seiner a-posteriori-Vorhersagefehlerquadrate F(k), eines Ruckwartslinearvorhersagekoef- 
fizientenvektors b(k) des Eingangssignals x(k) und der Summe seiner a-posteriori-Vorhersagefehlerquadua- 
drate B(k), und 

- Mittel (42, 43, 44) zum Ermitteln des Vorfilterkoeffizientenvektors g(k) durch eine Rekursionsfbrmel, die sich 
aus der folgenden ersten und zweiten Gleichung zusammensetzt: 



[V]- 



wobei f(k) einen Vorfilterableitungskoeffizientenvektor reprasentiert. 

10. Vorrichtung nach Anspruch 8 oder 9, bei der die Mittel zur Berechnung des inneren Produkts Korrelationsrechen- 
mittel zum Berechnen des Korrelationsvektors r^k) des Eingangssignals x(k) durch 

Vi( k > =r P-i( k - 1 )- x ( k - L )vi( k - L ) + x ( k )vi( k ) 

umfassen, wobei Xp.^k) den Vektor des Eingangssignals x(k) reprasentiert und ausgedruckt ist als: 

x M (k)=[x(k-1),x(k-2) x(k-p+1)] T 

11. Vorrichtung nach Anspruch 7 Oder 9, bei der das Vorfilterableitungskoeffizientenvektorkorrigiermittel Mittel sind, 
die, wenn das letzte Element des Vorfilterkoeffizientenvektors g(k-1) als durch g p (k-1) reprasentiert angenommen 
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wird, den Vorfilterkoeffizientenvektor g(k) unter Verwendung der folgenden Gleichung berechnen, die eine modffi- 
zierte Version der zweiten Gleichung ist: 



12. Vorrichtung nach Anspruch 7 oder 9, bei der, wenn zwei vorbestimmte nicht-negative Zahlen durch 6p(k) und ^(k- 
1) reprasentiert seien, die Vorfilterkoeffteientenvektorkorrigiermittel und die Vorfilterableitungskoeffizien- 
tenvektorkorrigiermittel den Vorfilterkoeffizientenvektor g(k) durch die folgenden Gleichungen berechnen, die modi- 
fizierte Versionen der ersten und der zweiten Gleichung sind: 



13. Vorrichtung nach Anspruch 7 oder 10, ferner umfassend Subtrahiermittel zur Bildung, als das Fehlersignal e(k), der 
Differenzy(k) - y (k) zwischen dem Ausgangssignal y(k) des unbekannten Systems (12) und dem geschatzten 
Signal y (k), welches das Ausgangssignal des Schatzsystems ist. 

Revendications 

1. Precede d'estimation adaptative de la fonction de transfert d'un systeme technique inconnu par utilisation d'un 
algorithme de projection, iedit precede comprenant : 

(i) ('application d'un signal physique d'entree x(k) audit systeme inconnu (12), k representant un temps discret ; 

(ii) la detection d'un signal physique de sortie y(k) provenant dudit systeme inconnu (12) ; 

(iii) I'estimation de la fonction de transfert en utilisant Iedit signal d'entree x(k) et Iedit signal de sortie y(k) et en 
generant un signal de sortie y (k) d'un systeme estime ayant ladite fonction de transfert estimee h (k), iedit sys- 
teme estime etant realise sous la forme d'un filtre numerique ayant un nombre de prises de L ; 

(iv) I'obtention, en tant que signal d'erreur e(k) = y(k) - y(k) , de la difference entre le signal de sortie y (k) dudit 
systeme estime et le signal de sortie y(k) dudit systeme inconnu (12) ; 

(v) le calcul d'un vecteur de correction 5h (k) pour la fonction de transfert estimee a partir de vecteurs de 

signaux d'entree x(k), x(k-1) x(k-p+1) et de signaux d'erreur e(k), e(k-1) e(k-p+1) , ou p est un entier 

egal ou superieur a 2 et x(k) est defini par x(k) = (x(k), x(k-1), x(k-L+1)] T ; et 

(vi) I'utilisation dudit vecteur de correction 8h (k) et d'une taille de pas de correction predeterminee u pour cor- 
riger de facon repetee ladite fonction de transfert estimee h (k) a chaque instant k conformement a I'equation 
suivante : 



afin que Iedit signal d'erreur e(k) tende vers zero ; 
caracterise en ce que I'etape (y) comprend les etapes : 

(v-1) de calcul d'un vecteur de coefficients de pre-filtre g(k) en resolvant le systeme d'equations lineaires 
a p inconnues suivant : 






h(k+1) = h(k)+u8h(k) 



R p (k)g(k) = e(k), 
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ou e(k) represente le vecteur dudit signal d'erreur e(k) et est defini par : 

e(k) = [e(k), (1-u)e(k-1) (l-^^efk-p+l)] 1 , et 

5 Rp(k) est une matrice de covariance dudit signal d'entree x(k) definie par : 

R p (k) = [x(k), x(k-1) x(k-p+1)] T [x(k), x(k-1) x(k-p+1)] ; et 

(v-2) d'utilisation dudit vecteur de coefficients de pre-filtre g(k) pour calculer ledit vecteur de correction 
w 8h (k) confbrmement A Pequation suivante : 

8h(k) = [x(k), x(k-1) x(k-p+1)]g(k), 

ou x(k) represente le vecteur dudit signal d'entree x(k) ; I'etape (v-1) comprenant les etapes : 
15 (v-1-a) de calcul d'un vecteur de coefficients de prediction lineaire vers I'avant a(k) dudit signal d'entree 

x(k), de la somme de ses carres de Perreur de prediction a posteriori F(k), d'un vecteur de coefficients de 
prediction lineaire vers Tamere b(k) dudit signal d'entree x(k) et de la somme de ses carres de I'erreur de 
prediction a posteriori B(k), et 

(v-1-b) posant qu'un vecteur de coefficients de derivation de pre-filtre est represente par f(k), d'obtention 
20 du vecteur de coefficients de pre-filtre g(k) par une formuie recursive constituee des premiere et seconde 

equations suivantes : 



[ 0 l a(k) T e(k) 

g(W - (i-4f(k-uJ + F ( k) a < k > 

fftk-Dl b(k-l) T e(k-l) 

L o J = *l k -U ~ B(FT) b(k 



35 2. Procede d'estimation adaptative de la fbnction de transfert d'un systeme technique inconnu par utilisation d'un 
algorithme de projection, ledit procede comprenant : 

(i) I'application d'un signal physique d'entree x(k) audit systeme inconnu (12), k representant un temps discret ; 

(ii) la detection d'un signal physique de sortie y(k) provenant dudit systeme inconnu (12) ; 

40 (iii) I'estimation de la fbnction de transfert en utilisant ledit signal d'entree x(k) et ledit signal de sortie y(k) et en 

generant un signal de sortie y (k) d'un systeme estime ayant ladite fonction de transfert estimee h (k), ledit sys- 
teme estime etant realise sous la forme d'un filtre numerique ayant un nombre de prises de L ; 
(iv) I'obtention, en tant que signal d'erreur e(k) = y(k)- y(k) , de la difference entre le signal de sortie y (k) dudit 
systeme estime et le signal de sortie y(k) dudit systeme inconnu (12) ; 

45 (v) la repetition du calcul de la fonction de transfert k chaque instant k afin que ledit signal d'erreur e(k) tende 

vers zero ; 

caracterise en ce que 

I'etape (y) comprend les etapes : 

so (v-1) de calcul d'un vecteur de coefficients'de pre-filtre g(k) en resolvant le systeme d'equations lineaires 

k p inconnues suivant : 

R p (k)g(k) = e(k), 

55 ou e(k) represente le vecteur dudit signal d'erreur e(k) et est defini par : 

e(k) = [e(k). (1^)e(k-1) (l-nJ^etk-p+l)] 1 , et 
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R p (k) est une matrice de covariance dudit signal d'entree x(k) d6finie par : 
R p (k) = Wk), x(k-1) x(k-p+1)] T [x(k), x(k-1) x(k-p+1)] ; 

(v-2) de (issage des coefficients de pre-filtre grfk) qui sont des 6!6ments d'un vecteur de coefficients de 
pre-filtre g(k) confbrmement aux Equations suivantes : 

Si(k) = s M (k-1)+ug rfk) pour 2<;i<sp 
= U9i(k) pouri=1 

afin d'obtenir un coefficient de iissage s,(k) ; 

(v-3) d'obtention d'une fonction de transfert estimee approximative z(k+1) en tant qu'approximation de 
ladite fonction de transfert estimee en utilisant te coefficient de Iissage s p (k) sur la base de liquation sui- 
vante : 

2(k+1) = z(k) + s p (k)x(k-p+1); 

(v-4) de calcul de la convolution x(k) T z(k) de ladite fonction de transfert estimee approximative z(k) avec 
un vecteur de signal d*entr6e x(k) defini par x(k) = [x(k), x(k-1) x(k-L+1)] T ; 

(v-5) de calcul du produit scalaire s p . 1 (k-1) T r ^(k) , en fixant le vecteur dudit coefficient de Iissage Sj(k) 
et un vecteur de correlation r^k) dudit signal d'entree respectivement d : 

s p . 1 (k-1) = [s 1 (k-1),s 2 (k.1) s M (k.1)] T , 

et 

Vi(k) = Wk) T x(k-1), x(k) T x(k-2) x(k) T x(k-p+1)] T ; et 

(v-6) d'obtention de la somme dudit r6sultat de convolution x(k) T z(k) et dudit produit scalaire en tant que 
ledit signal de sortie estim£ y (k). 

Procede selon la revendication 2, dans lequel retape (v-1) comprend : 

(a) le calcul d'un vecteur de coefficients de prediction iineaire vers I'avant a(k) dudit signal d'entree x(k), de la 
somme de ses carr6s de I'erreur de prediction a posteriori F(k), d'un vecteur de coefficients de prediction 
Iineaire vers I'arriere b(k) dudit signal d'entree x(k) et de la somme de ses carres de I'erreur de prediction a 
posteriori B(k), et 

(b) posant que le vecteur de coefficients de derivation de pr6-filtre est represente par f(k), I'obtention dudit vec- 
teur de coefficients de pre-filtre g(k) par une formule recursive constitu6e des premiere et seconde equations 
suivantes : 



g(k) = (l 

["V 1 '] - 



Precede selon la revendication 1 ou 3, dans lequel, posant que le dernier element dudit vecteur de coefficients de 
pr6-filtre g(k-1) est represente par g p (k-1). ledit vecteur de coefficients de pre-filtre g(k) est calcule en utilisant 
I'equation suivante qui est une version modifiee de ladite seconde equation : 



Jo] a(k) T e(k) 
"^Lf (k-l)J + F(k) *< k > 

b(k-l) T e(k-l) 

: 9 (k - 1J 5TFT) b <*-D 
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= g(k-l) - g P (k-l)b(k-l) 



Precede selon la revendication 1 ou 3, dans lequel, posant deux nombres non negatifs predetermines comme 
etant represents par 5p(k) et ^(k-l), (edit vecteur de coefficients de pre-filtre g(k) est calcule par les equations 
suivantes qui sont des versions modifiees desdites premiere et seconde equations : 

g<k) = (i-^,[ f ( £ J + a(k, 

r f(k_1 )l b(k-l) T e(k-l) 

L 0 J - " B(k-l)^(k-l) *>l*-U 

Procede seion la revendication 2 ou 3, dans lequel I'etape (v-5) comprend le calcul dudit vecteur de correlation Tp_ 
^k) dudit signal d'entree x(k) conformement a : 

ViW = r p . 1 (k-1)-x(k.L)x p . 1 (k-L)+x(k)x p . l (k) 

ou x p . 1 (k) represente le vecteur dudit signal d'entree x(k) exprime par : 

x p . 1 (k) = [x(k-1),x(k-2) x(k-p + 1)] T 

Dispositif d'estimation adaptative de la fbnction de transfert d'un systeme technique inconnu par utilisation d'un 
algorithme de projection, led it dispositif comprenant : 

un moyen pour appliquer un signal physique d'entree x(k) audit systeme inconnu (12), k representant un temps 
discret ; 

un moyen pour detecter un signal physique de sortie y(k) provenant dudit systeme inconnu (12) ; 
un moyen (11) pour estimer la fbnction de transfert en utilisant ledit signal d'entree x(k) et ledit signal de sortie 
y(k) et en gen era nt un signal de sortie y (k) d'un systeme estime (23) ayant ladite fonction de transfert estimee 
h (k), ledit systeme estime etant realise sous la forme d'un filtre numerique ayant un nombre de prises de L ; 
un moyen (24) pour obtenir, en tant que signal d'erreur e(k) = y(k) - y(k) , la difference entre le signal de sortie 
y (k) dudit systeme estime et le signal de sortie y(k) dudit systeme inconnu (12) ; 

un moyen (21) pour calculer un vecteur de correction 5h (k) pour la fonction de transfert estimee a partir de 

vecteurs de signaux d'entree x(k), x(k-1) x(k-p+1) et de signaux d'erreur e(k), e(k-1) e(k- p+1) , ou p 

est un entier egal ou superieur a 2, et x(k) est defini par x(k) = [x(k), x(k-1) x(k) x(k-L+1)] T ; et 

un moyen (22) pour corriger de facon repetee, a chaque instant k, en utilisant ledit vecteur de correction 5h (k) 
et une taille de pas de correction predetermined u., ladite fonction de transfert estimee h (k) conformement a 
I'equation suivante : 

h(k+1) = h(k)+u8 h(k) 
afin que ledit signal d'erreur e(k) tende vers zero ; 

caracterise en ce que ledit moyen (21) de calcul du vecteur de correction 5h (k) comprend : 

un moyen (31) pour calculer un vecteur de coefficients de pre-filtre g(k) en resolvant le systeme d'equa- 
tions lineaires.a p inconnues suivant : 

R p (k)g(k) = e(k), 
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ou e(k) represente le vecteur dudit signal d'erreur e(k) et est defini par : 

e(k) = [e(k), (1-|i)e(k-1) (1-u) M e(k-P+1)] T . et 

5 Rp(k) est une matrice de covariance dudit signal d'entree x(k) definie par : 

R p (k) = [x(k), x(k-1) x(k-p+1)] T [x(k), x(k-1) x(k-p+1)] ; et 

un moyen (32) pour calcuier, en utilisant ledrt vecteur de coefficients de pre-filtre g(k), ledit vecteur de cor- 
10 rection 5h (k) conformement a r equation suivante : 

8h(k) = [x(k), x(k-1) x(k-p+1)]g(k) 

ou x(k) represente le vecteur dudit signal d'entree x(k) ; et 

ledit moyen (31) de calcui dudit vecteur de coefficients de pre-filtre g(k) comprend 



15 



un moyen (41) pour calcuier un vecteur de coefficients de prediction lineaire vers I'avant a(k) dudit signal 
d'entree x(k), la somme de ses carres de I'erreur de prediction a posteriori F(k), un vecteur de coefficients 
20 de prediction lineaire vers I'arriere b(k) dudit signal d'entree x(k) et la somme de ses carres de I'erreur de 

prediction a posteriori B(k), et 

un moyen (42, 43, 44) pour obtenir ledit vecteur de coefficients de pre-filtre g(k) par une formule recursive 
constitute des premiere et seconde equations suivantes : 



ou f(k) represente un vecteur de coefficients de derivation de pre-filtre. 

35 8. Dispositif d'estimation adaptative de la fonction de transfert d'un systeme technique inconnu par utilisation d'un 
algorithme de projection, ledit dispositif comprenant : 

un moyen pour appliquer un signal physique d'entree x(k) audit systeme inconnu ( 1 2), k representant un temps 
discret ; 

un moyen pour detecter un signal physique de sortie y(k) provenant dudit systeme inconnu (12) ; 
un moyen (11) pour estimer la fonction de transfert en utilisant (edit signal d'entree x(k) et ledit signal de sortie 
y(k) et en generant un signal de sortie y (k) d'un systeme estime (23) ayant ladite fonction de transfert estimee 
h (k), ledit systeme estime etant realise sous la forme d'un filtre numerique ayant un nombre de prises de L ; 
un moyen (24) pour obtenir, en tant que signal d'erreur e(k) = y(k) - y(k) , la difference entre le signal de sortie 
y (k) dudit systeme estime et le signal de sortie y(k) dudit systeme inconnu (12) ; et 

un moyen (20) pour calcuier la fonction de transfert estimee a chaque instant k afin que ledit signal d'erreur 
e(k) tende vers zero ; 

caracterise en ce que ledit moyen (20) de calcui de ladite fonction de transfert estimee comprend : 

so - un moyen (31) pour calcuier un vecteur de coefficients de pre-filtre g(k) en resolvant le systeme d'equa- 

tions lineaires a p inconnues suivant : 

R p (k)g(k) = e(k), 

55 ou e(k) represente le vecteur dudit signal d'erreur e(k) et est defini par : 

e(k) = [e(k), (1-u)e(k-1) (l-uj^etk-p+l)] 1 , et 



40 
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Rp(k) est une matrice de covariance et ledit signal d'entree x(k) est d6fini par : 
R p (k) = Wk), x(k-1) x(k-p+1)] T [x(k), x(k-1) x(k-p+1)] ; et 

un moyen (51) pour lisser les coefficients de pre-filtre gtfk). qui sont des elements dudlt vecteur de coeffi- 
cients de pre-filtre g(k), conformement aux equations suivantes : 

s,(k) = s M (k-1)+|ig i (k) pour2^p 
= MiW pour 1=1 



afin d'obtenir un coefficient de lissage Sj(k) ; et 

- un moyen (55) pour obtenir une fonction de transfert estimee approximative z(k+1) en tant qu'approxima- 
tion de ladite fonction de transfert estimee en utilisant le coefficient de lissage s p (k) sur la base de liqua- 
tion suivante : 

z(k+1) = z(k) + s p (k)x(k-p+1);et ' 
ledit systeme estime (23) comprend : 

- un moyen (54) pour calculer la convolution x(k) T z(k) de ladite fonction de transfert estimee approximative 
z(k) avec un vecteur de signal d'entree x(k) defini par x(k) = [x(k), x(k-1), .... x(k-L+1)] 

- un moyen (53) pour calculer le produit scalaire s ^ (k-1 ) V ^ (k) , en fixant respectivement le vecteur dudit 
coefficient de lissage Sj(k) et un vecteur de correlation r^k) dudit signal d'entree, d : 

Sp.^k-1) = [s ^k-1), s 2 (k-1) s^fk-l)] 7 , 

et ' 

r M (k) = [x(k) T x(k-1), x(k) T x(k-2) x(k) T x(k-p+1)] T ; et 

- un moyen (57) pour obtenir la somme dudit resultat de convolution x(k) T z(k) et dudit produit scalaire en 
tant que signal de sortie estime y (k). 

9. Dispositif selon la revendication 8, dans lequel ledit moyen (31) pour resoudre ledit systeme d'equations lineaires 
a p inconnues comprend : 

un moyen (41) pour calculer un vecteur de coefficients de prediction lineaire vers I'avant a(k) dudit signal 
d'entree x(k). la somme de ses carres de I'erreur de prediction a posteriori F(k), un vecteur de coefficients de 
prediction lineaire vers I'arriere b(k) dudit signal d'entree x(k) et la somme de ses carres de I'erreur de predic- 
tion a posteriori B(k), et 

un moyen (42, 43. 44) pour obtenir ledit vecteur de coefficients de pre-filtre g(k) par une formule recursive 
constituee des premiere et seconde equations suivantes : 

ou f(k) represente un vecteur de coefficients de derivation de pre-filtre. 

10. Dispositif selon la revendication 8 ou 9, dans lequel ledit moyen de calcul de produit scalaire comporte un moyen 
de calcul de correlation pour calculer ledit vecteur de correlation r^k) dudit signal d'entree x(k) conformement & : 

'p-iM = r p . l (k-1)-x(k-L)x p . l (k-L)+x(k)x p . 1 (k) 
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ou Xp.i(k) represents le vecteur dudit signal d'entr6e x(k) exprim6 par : 



x M (k) = [x(k-1),x(k-2) x(k-p+1)] T 



11. Dispositif selon la revendication 7 ou 9, dans lequel ledit moyen de correction de vecteur de coefficients de deriva- 
tion de pr£*fittre est un moyen qui, posant que le dernier element dudit vecteur de coefficients de pr6-filtre g(k-1) 
est represents par g p (k-1), calcule ledit vecteur de coefficients de pre-filtre g(k) conform 6ment & liquation sun 
vante, qui est une version modtfi6e de ladite seconds equation : 



12. Dispositif selon la revendication 7 ou 9, dans lequel, posant que des nombres non negattfs predetermines sont res- 
pectivement representes par 6p(k) et 8b(k-1), ledit moyen de correction de vecteur de coefficients de pr6-filtre et 
ledit moyen de correction de vecteur de coefficients de derivation de pr6-filtre calculent ledit vecteur de coefficients 
de pre-filtre g(k) respectivement & partir des equations suivantes, qui sont des versions modifi6es desdites pre- 
miere et seconde equations : 



13. Dispositif selon la revendication 7 ou 10, comprenant en outre un moyen de soustraction pour fournir, en tant que 
ledit signal d'erreur e(k), la difference y(k)- y (k) entre ledit signal de sortie y(k) dudit systeme inconnu (12) et ledit 
signal estime y (k) qui est la sortie dudit systeme estime. 
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